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■ Abstract. Gravitational waves contain tail effects that are due to the backscattering 

of linear waves in the curved space-time geometry around the source. The knowledge 
as well as the accuracy of the two-body inspiraling post-Newtonian (PN) dynamics and 
of the gravitational-wave signal has been recently improved, notably by computing the 
spin-orbit (SO) terms induced by tail effects in the gravitational-wave energy flux at the 
$H . 3PN order. Here we sketch this derivation, which yields the phasing formula including 

SO tail effects through the same 3PN order. Those results can be employed to improve 
the accuracy of analytical templates aimed at describing the whole process of inspiral, 
merger, and ringdown. 

> : 

1. Introduction 

d. 

By 2016 - 2018, the ground-based gravitational-wave detectors Virgo and LIGO will 
be upgraded to such a sensitivity that event rates for coalescing binary systems will 
increase by approximately a factor one thousand, making likely the first detection of 
gravitational waves from those systems. The search for gravitational waves from coa- 
lescing binary systems and the extraction of source parameters require a rather accurate 
knowledge of the waveform of the incoming signal. The post-Newtonian (PN) ex- 
pansion is the most powerful approximation scheme in analytical relativity capable of 
describing the two-body dynamics and the gravitational-wave emission of inspiraling 
compact binary systems (Blanchet 2006). 

The presence of spin effects adds substantial complexity to the gravitational wave- 
forms, making it indispensable to include them into search templates. We have recently 
improved (Blanchet et al. 2011) the knowledge and accuracy of the two-body inspi- 
raling dynamics and gravitational-wave signal by computing the spin-orbit (SO) terms 
induced by tail contributions due to the back-scattering of linear waves in the curved 
space-time geometry around the source. Here we shall summarize this work, which is 
the continuation of Faye et al. (2006) and Blanchet et al. (2006) where we obtained the 
2.5PN SO contributions in the equations of motion and gravitational-wave energy flux. 

After briefly reviewing the post-Newtonian multipole moment formalism and dis- 
cussing relevant properties of tails, we describe how spin effects are included to our 
scheme and derive the binary's evolution equations when black holes carry spins. Next, 
we investigate the time evolution of the moving triad and solve the precessing dynamics 
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at the relevant PN order. We then compute the tails, which depend on the recent past 
history of the source, restricting ourselves to quasi-circular adiabatic inspiral. At last, 
we derive in the 3PN SO tail effects in the energy flux and in the gravitational phasing. 

2. Wave generation formalism 

The gravitational waveform hJJ, generated by an isolated source described by a stress- 

I J 

energy tensor T MV with compact support, and propagating in the asymptotic regions 
of the source, is the transverse trace-free (TT) projection of the metric deviation at 
the leading-order l/R in the distance to the source. It is parametrized by symmet- 
ric trace-free (STF) mass-type moments Ul and current-type ones Vl, referred to as 
radiative moments, which constitute observable quantities at infinity from the source. 
The general expression of the TT waveform, in a suitable radiative coordinate system 
= (c T, X), reads (Thorne 1980), when neglecting terms 0(1 /R 2 ) with R = |X| 



The moments Ul and Vl are functions of the retarded time Tr = T - R/c. The integer 
I refers to the multipolar order and N = X/R is the unit vector pointing from the 
source to the far away detector. The tensor Pj^ denotes the transverse-traceless (TT) 

projector PnuPi)j - PijPul% where Pjj - djj - N,Nj is the projector orthogonal to 
N. The quantity is the Levi-Civita symbol such that e 12 3 = 1. Round parentheses 
indicate symmetrization. After plugging Eq. (1) into the standard expression for the 
gravitational-wave energy flux, we get 



where the superscript (1) stands for the first time derivative (Thorne 1980). 

In the multipolar-post-Minkowskian formalism, developed by Blanchet & Damour 
(1986, 1988, 1992), the radiative moments are linked to six sets of multipole moments 
characterizing the source, collectively called the "source" moments and denoted I L , 
Jl, Wl, Xl, Yl, Zl. The relation between Ul, Vl and the former quantities encode 
all the non-linearities in the wave propagation between the source and the detector. 
After being re-expanded in a PN way, they are seen to contain the contribution of the 
so-called gravitational-wave tails, due to backscattering of linear waves onto the space- 
time curvature associated with the total mass of the source itself. The expression for 
Ul(Tr) at the 1.5PN order, where the tail effects first appear, is 



Here t' is a freely specifiable time scale (Blanchet 1995). A similar equation links 
Vl(Tr) to J L . For the present application, Wl, X l , Y l and Z L , associated to a possible 
gauge transformation performed at linear order, play no role. It will be sufficient to 
consider II at 1.5PN and Jl at 0.5PN order. The 1.5PN mass moments are given by 




(1) 




(2) 




(3) 
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an integral extending over the mass density cr = (T 00 + T")/c 2 and the current density 
cTj = T°'/c of the matter source (Blanchet et al. 2006): 



r , 3 i . , l2 (2) 4(2^+1) , 

= d x x L cr + — xl x cr ; xil cr. + (9 — 

J L 2c 2 (2^ + 3) c 2 (^ + 1)(2£ + 3) ' J V 4 ' 



(4) 



The 0.5PN part of Jl is obtained by substituting £ a b(i c XL-\)a o~b to %l cr in the first term. 

To control the past behavior of the tail integral in Eq. (3), we assume that, at 
early times, the source was formed from a bunch of freely falling particles moving on 
some hyperbolic-like orbits, and forming at a later time a gravitationally bound sys- 
tem by emission of gravitational radiation. This ensures that the integral UiSTr) = 

At l L t+2 \t) ln[(7# - 0/2Tq] is convergent. To compute it, it is convenient to perform 

the Fourier decomposition Iiif) = dQ./(2n) e~ lQt , and commute the tail inte- 
gral with the Fourier one. We obtain a closed-form expression by resorting to standard 
mathematical formulae for the integrals of e~ lSlt In t. The result reads (with s(Q) = || 
and je being the Euler constant) 



— I L (n) e~ iaT « [-s(Q) + i(ln(2|Q|r ) + y E ) 



(5) 



3. Applications to spinning binaries 
3.1. Spin vectors for point-like objects 

Following our previous works (Faye et al. 2006; Blanchet et al. 2006), we base our 
calculations on the model of point-particles with spins (Mathisson 1937; Papapetrou 
1951; Tulczyjew 1959; Bailey & Israel 1980). The stress-energy tensor T /iy of a system 
of spinning particles is the sum of a monopolar piece, made of Dirac delta-functions, 
plus the dipolar or spin piece, made of gradients of delta-functions: 

where £ (4) is the four-dimensional Dirac function, x^ is the field point, is the world- 
line of particle A, u A - dy^/(cdTA) is the four- velocity such that g^ v u^u A - -1 (with 
gfiv = Sfiv(yA) denoting the metric at the particle's location), p^ A is the linear momentum 
of the particle, and denotes its antisymmetric spin angular momentum. Our spin 
variables has been rescaled by a factor c (S^ v = cS^ e ) so as to have a non-zero New- 
tonian limit for compact objects. Imposing that the antisymmetric part of T MV should 
vanish yields the covariant equations of evolution for the spin, while the covariant equa- 
tions of motion follow from the conservation relation V y r /JV = (Papapetrou 1951). 

In order to fix unphysical degrees of freedom associated with an arbitrariness in 
the definition of S^ v in the case of point particles (due to the freedom in the choice 
of the location of the center-of-mass worldline within the extended bodies), we adopt 
the covariant supplementary spin condition S^Py = (Tulczyjew 1959). We restrict 
the computation to linear terms in the spins, neglecting 0(S A ). Then the mass = 
(.-p^pft) 112 is conserved and S^ A is parallel transported along the worldline of body A. 



cPj" ,,v 
^ A U A 



d^\x - y A ) 



(6) 
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As is standard, we introduce a spin vector variable S' A with constant magnitude (see 
e.g., Kidder 1995), for instance by projecting onto some orthonormal basis of the 
A-particle rest space and taking the Hodge dual with s,-^. 

The post-Newtonian expressions of the densities cr, cr,- entering the source mo- 
ments (4), as well as the stress density o",y = T'-i, are obtained iteratively from the 
truncated components of the stress-energy tensor (6). They are used in turn to compute 
the metric at the next iteration step of the perturbative scheme. 

Inserting the former matter densities into Eq. (4) (or its counterpart for V£), we 
see that spins arise at 1.5PN order in li and 0.5PN order in J^. By virtue of Eq. (2), 
leading SO terms in the gravitational- wave flux are thus seen to be of half-integer PN 
order. However, due to the 1 /c 3 factor in front of the tail integrals, as shown in Eq. (3), 
terms of integer orders may also appear: (i) when the source moments are contained 
in the tail integrands, or (ii) when they are multiplied by non-spin tail integrals. The 
leading 3PN terms coming from C/,y and Vy are precisely those we want to compute. 

3.2. 1.5PN dynamics with spin-orbit effects 

To reduce the accelerations generated by the time differentiations of /,y and Jjj and, 
most importantly, to find the time dependence of if) and jf) required for calculating 
integrals such as the one in Eq (3), we shall rely on the 1.5PN dynamics. The corre- 
sponding PN equations of evolution at this order are derived by inserting the 1.5PN 
metric into the co variant equations and discarding (9(1 /c 4 ) remainders. Going to the 
center-of-mass frame, we denote by x = yi - y2 the relative position of the particles 
(and r = |x|, v - dx/dt). Following Kidder (1995) we introduce an orthonormal moving 
triad {n, A, {} defined by n = x/r, t = L N /|L N |, where L N = mv x x v denotes the Newto- 
nian orbital angular momentum (with m - m\ + m2 and v - m\m2lm 2 ), and A - i X n. 
For general orbits, the time derivatives of those vectors may be written in the form 



dn y dA dt 

— - CO A , — = -OJ11- CO prec I , — = CO pmc A . (7) 

dt dt dt 



The angular frequency co at which the separation n rotates in the instantaneous orbital 
plane is the so-called orbital frequency, while the third equation (7) defines the preces- 
sional frequency a» prec , which gives thus the variation of I in the direction of A. 

Instead of writing the equations of motion in the form of a link between d\/dt and 
the kinematic variables, we give the expressions of co, co piec . Focusing from now on to 
the case of quasi-circular motion, for which r is constant apart from the 2.5PN secular 
radiation damping, we have (see e.g., Faye et al. 2006) 

co 2 = -^{l + y (-3 + v) + y 3/2 (-5^ - 36cr e )j + 0(-) , (8) 
wp re c - -oo y 3/2 (ls n + 36cr n ) + tf(-) , (9) 

where the dimensionless spin variables are defined by s = (Si + S2) /(Gm 2 ) and cr - 
($2/11*2 ~ Si/mi)/(Gm), whereas 8 = X\ - X2, S( = s • t and s n = s • n; y = Gm/(rc 2 ) is 
the PN parameter. The spin vectors Sa satisfy the usual-looking precession equations 



= £Ia x Sa , with il\ - coy 

dt 



3 v 3 

4 + 2-4 5 



do) 



3PN spin-orbit effect in the gravitational radiation ffux of compact binaries 



5 



showing that the spins, at the 1.5PN order, precess around the direction of (, and at the 
quasi-constant rate Qa = \&aI To obtain Q.2 we simply have to change 6 into -8. 



4. Computation of the tails 

The only time dependence in the source moments composing the tail integrals occurs 
through the triad vectors and the spins at the 1.5PN order. We adopt the following 
strategy to integrate Eqs. (7): (i) we express the three moving triad vectors in terms of 
the Euler angles that link the latter triad to some orthonormal fixed triad whose third 
vectors points to the direction of the total angular momentum at some time to; (ii) we 
express the appropriate brick combinations of those angles as a function of S„, S^, Sf 

and <f> = J dtoj(t); (iii) we solve for the precession equations (10) in the triad basis and 
insert the results into the formulae we got for n, A, ( in the previous step. 

Although the binary continuous spectrum of frequencies of li{t) should contain 
all orbital frequencies at any epoch in the past, say oj(t) with t < Tr, differing from the 
current orbital frequency o)(Tr) due to gravitational radiation damping, we can actually 
compute the tail integral by considering only the current frequency o)(Tr), modulo 
small error terms of negligible order 0(ln c/c 5 ). Similarly, we can regard the precession 
frequencies J2a(0 of the two spins as constant in the calculation. Thus, we may replace 
<p(t) by co(t - to) + <p(to) and Qa(0 by CIa(Tr) in the tail integrand. This takes the form 
of a sum of complex exponentials of the type g 1& W, with o) npq = no) + pQ.\ + qQ.2, 
times certain coefficients A^ q . It is then straightforward to take the Fourier transform 
of the source moments, e.g. Iiif), and apply formula (5), which leads to 



Ul(T r ) = ^ iA n L P H- 

n,p,q 



■10). 



npq) & 



^s(o) npq ) + i(ln(2|w„ M |To) + te) 



(11) 



5. Energy flux and orbital phasing 

Inserting the explicit expressions for U,j and Vy into the energy flux (2) yields the 
net SO tail contribution at 3PN order, in terms of the gauge-invariant PN parameter 
x = (Gma>/c 3 ) 2/3 and the spins. The SO terms have to be included in either the in- 
stantaneous moment in front of the tail integrals, or in the tail integrals themselves. 
This gives several "direct" SO contributions coming from tails at relative 1.5PN order 
(for the mass quadrupole tail) or 0.5PN order (for the current quadrupole tail) which 
are then added together. In addition there is the crucial contribution due to the reduc- 
tion to circular orbits of the standard (non-spin) tail integral at 1.5PN order, for which 
the relation between the orbital separation r and the orbital frequency o) [as given by 
the inverse of Eq. (8)] provides a supplementary SO term at relative 1.5PN order, which 
contributes in fine at the same 3PN level as the "direct" SO tail terms. Finally we obtain 

6r = ^ C l x W[-l6nse- 3 -^6cr { }, (12) 

Using an energy balance argument, we finally equate the averaged evolution of the 
binding energy of the binary reduced for quasi-circular orbits dE(x, S(, 07)/d? (see e.g., 
Blanchet et al. 2006) to -T(x, st, ere), where T is the sum of the 2PN SO flux obtained 
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by Kidder (1995), the 3PN non-spin contribution investigated by Blanchet et al. (2002, 
2004) and the correction (12) we just calculated. The time derivative applied to E 
acts only to the frequency variable since, as we have checked, S{ and 07 are secularly 
conserved (neglecting SS contributions). We rind the following secular variation of 
frequency, denoted 6j for simplicity, at the 3PN order including SO effects: 

^ 96 ,,,r / 743 11 v 47 25 v 

2 /34103 13661 59 2 x 

+X (l8144 + 20l6" y+ T8 V ) 

,„/ 4159 5861 809 r 189 1001 281 n x 

+j — — — 7T — 5/ - — -<5cr> + v — — n + — —st + -t-SctA) 

V 672 144 84 L g 12 8 " 

,/ 16447322263 16 , 1712 856 188/r 151tt 

+X ( 139708800 + T" -l05 7E -W5 lnil6x) -— Si -— 6(re 
r 56198689 451 2n 541 , 5605 ,^ 

+y [-^1772T + 18"" 1 + 896 V - 2592 V )) • (13) 

Integrating this, we get the "carrier" phase 0(?) = fiodt. The total phase O = - a + 
0(\/c 4 ), where a is the precessional phase, can be computed numerically or manually. 
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